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O I Abstract 
^ . 

■ Using a recently developed perturbation formalism based on curvature 

I quantities, we investigate the linear stability of black holes and solitons 
with Yang-Mills hair and a negative cosmological constant. We show 

^ I that those solutions which have no linear instabilities under odd- and 

bJQi even-parity spherically symmetric perturbations remain stable under odd- 

^ I parity, linear, non-spherically symmetric perturbations. 
. 

^ ! 1 Introduction 

The discovery of solitonic [Q] and black hole ||2| solutions to 5u(2) Einstein- 
Yang- Mihs (EYM) theory sparked considerable study of the properties of non- 
Abelian gauge theories coupled to gravity (see |^ for a comprehensive review 
of the subject and a full list of references). Many examples of both globahy 
regular and black hole solutions have now been discovered, in a wide range of 
theories involving various non-Abehan gauge fields, both with and without a 
(positive or negative) cosmological constant. 

Of this plethora of examples, many (but not all) are unstable. In particular, 
there is a general result that solitons and black holes solutions of EYM theory 
with a compact gauge group in asymptotically flat space must be topologically 
unstable Q. This instability is analogous to the instability of the flat-space 
Yang-Mills-Higgs sphaleron |p. In addition, su(2) EYM solitons and black 
holes in the presence of a positive cosmological constant are unstable It 
was therefore surprising to discover that this result does not extend to solutions 
in anti-de Sitter (adS) space, when there is a negative cosmological constant. 
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Both soliton Q and black hole |^ solutions have been found which are linearly 
stable with respect to spherically symmetric perturbations. It is the purpose 
of the present article to examine the linear stability of these solutions with 
respect to non-spherical perturbations. 

We concentrate here on the odd-parity (or sphaleronic) sector. For spheri- 
cally symmetric perturbations, the behaviour in this sector is "topological" in 
the sense that it does not depend on the detailed structure of the solutions, but 
only on global properties and boundary conditions. Further, it is in this sector 
(which, for spherically symmetric perturbations, involves the perturbations of 
the gauge field only, and not the metric perturbations) in which the analogy 
with the flat space sphaleron is most pertinent. Analysis of this sector for soli- 
tons and black holes in asymptotically flat space has shown that all the modes 
of instability are contained within the spherically symmetric perturbations |^]. 
Again, this is the same as for the flat space sphalerons, where the only modes of 
instability are spherically symmetric [0]. In this article we extend this result 
to solitons and black holes in asymptotically anti-de Sitter space, and show 
analytically that those solutions which are stable under spherically symmetric 
perturbations remain stable under non-spherically symmetric perturbations in 
the odd-parity sector. We shall return to the question of the stability under 
non-spherically symmetric perturbations in the even-parity sector in a future 
publication. 

In order to discuss the stability, we use a recently developed perturbation 
formalism, which is based on curvature quantities. The main advantage of this 
formalism is that it allows us to cast the pulsation equations, governing linear 
fluctuations on a static background, into the form of a gauge invariant wave 
equation even when complicated matter fields are coupled to the metric. It has 
been shown in that for a static and purely magnetic solution of the EYM 
equations (with arbitrary compact gauge group), the pulsation equations admit 
the form of a symmetric wave equation for the linearized extrinsic curvature 
and the electric field. It is precisely this symmetric form of the pulsation 
equations which will be crucial in this article to show the stability of the above 
mentioned solitons and black holes by analytic means. 

This work is organized as follows. In section ^ we remember some impor- 
tant results about the hairy black hole solutions with a negative cosmological 
constant, recently found in |^], and the corresponding solitonic solutions 0. 
In section ^, we briefly review the curvature-based formalism of perturbation 
theory for a static background, and also generalize to the case where a cosmo- 
logical constant is present. The harmonic decomposition and the decoupling of 
the constraint and dynamical variables are performed in section ^, where the 
absence of exponentially growing modes is also shown. Technical details and 
the initial value formulation in terms of gauge-invariant quantities are discussed 
in Appendix ^ and respectively. 

The metric signature is (—,+,+,+) throughout, and we use the standard 
notations 2iO(^ab) = ^ab + ^ba and 2u;[ab] = ^ab — ^ba for symmetrizing and 
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antisymmetrizing, respectively. Throughout the paper, greek letters denote 
spacetime indices taking values in (0,1,2,3), while roman letters will denote 
spatial indices taking values (1,2,3). 



2 Solitons and hairy black holes with a negative cos- 
mological constant 

Hairy black holes in 5u(2) Einstein- Yang-Mills theory with a negative cosmo- 
logical constant were first found in j^, and subsequently in [^, where the 
corresponding solitonic solutions are also discussed. The equilibrium metric is 
spherically symmetric 

ds^ = -N{r)S^{r) dt^ + N-\r) dr^ + {dO^ + sin^ 9 dcj?) , 

and the gauge field potential has the spherically symmetric form 

A = {1- w{r)) [-T^dO + Te sin 6 . 

Here the su(2) generators Trfi^^f, are given in terms of the usual Pauli matrices 
C7j by Tr = Cj. ■ g_/2i, etc. Writing N{r) = 1 — 2m{r)/r — Ar^/3, where A is the 
(negative) cosmological constant, the field equations take the form: 



rrir = G 



Sr _ 2Gw^ 

'J ~ 



w 
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= Nr^Wrr + [2m G- — U;^ + (1 - w^)w, (1) 



where G is Newton's constant, and we have set the gauge coupling constant 
equal to \/47r for convenience. 

We are interested in solutions which approach anti-de Sitter (adS) space at 
infinity, so that the asymptotic behaviour of the field functions is: 

m{r) = M + — + 0(r~'^), 
r 

w[r) = Woo-\ \-0[r ), 



r 



S{r) = l+0{r 



Here we already observe one difference between the configurations in asymp- 
totically flat, and asymptotically adS space. For asymptotically flat space 
solutions, either Woo = 0, in which case the Reissner-Nordstrom (RN) solution 
follows, or Woo = if; so that there is no magnetic charge at inflnity [|l2[. How- 
ever, in adS, the boundary conditions place no restriction on the value of Woo, 
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so in general, even non-Abelian solutions will be globally magnetically charged. 
For black hole solutions having a regular event horizon at r = r/j, all the field 
variables have regular Taylor expansions near the event horizon, for example 

w{r) = w{rh) + Wr{rh){r - r/j) + 0{r - ruf . 

However, there are just two independent parameters, S{rh) and w{rh) since 
= at the event horizon, which gives 

^[rh) = y - — ■ 

In order for the event horizon to be regular, we shall also require that Nr{rh) > 
0, which implies that 

F,^l-Ar^-GM!:4^>0. 

From (1), one has 

_ {wjrhf - l)wirh) 
^"-^ '^ ~ r,F, 

There are also globally regular (solitonic) solutions, for which the behaviour 
near the origin is: 

m(r) = 2G6V3 + 0(r^), 
w{r) = l-6r^ + 0(r^), 
S{r) = S{0)[l + 4Gb^r^ + 0{r^)] . 

Here the independent parameters are b and S'(O). 

The simplest solutions to the field equations (|^) are the Schwarzschild-adS 
solution, 

w = 1, 5 = 1, m = const. 

and the RN-adS solution 

Q 

w = 0, S = 1, m = const. . 

2r 

In both cases, the YM field is effectively Abelian. It is however interesting to 
study the stability of these solutions with respect to non-Abelian perturbations. 

The solutions of greatest interest in this article are effectively non-Abelian 
solutions, for which the gauge function w has no zeros, since these solutions 
were shown in ||8| to be linearly stable to both even and odd-parity spherical 
perturbations. These solutions are unique to anti-de Sitter space, as w must 
have at least one zero if the cosmological constant is positive or zero 14] 



In Q it is proved that for any value of the gauge field at the event horizon, 
w{rh) 7^ 0, for all sufficiently large |A| there is a black hole solution in which w 
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has no zeros. Similar behaviour is found numerically for the solitonic solutions 
0]. For spherical perturbations of these equilibrium configurations, it is proved 
analytically that all the solutions in which w has no zeros are stable in the odd- 
parity sector 1^]. The even-parity sector is more complicated, but stability can 
be proven for sufficiently large |A|. 

In section ^ we shall require various properties of these soliton and black 
hole solutions, particularly when |A| ^ 1. We now briefly review the relevant 
results, and refer the reader to |8| for further details and proofs. We are 
concerned with black hole and soliton solutions in which the gauge field function 
w has no zeros, which exist for sufficiently large |A|. As |A| oo, the function 
w approaches a constant value, given by w{rh) for black hole solutions, and 
w = 1 for solitons. This means that Wr{r) tends to zero for all r as |A| ^ oo. 
In 1^, it was shown that in fact 

Wr{r) ~ o(|A|^^) as |A| oo 



for the black hole solutions, and that proof is easily extended to show that the 
same is also true for the solitons. This result will be useful in section |^. 

The other result we shall require in section |^ is that for these solutions in 
which w has no zeros, the gauge function w cannot be equal to ±1 anywhere, 
with the exception of the origin for the solitons (where regularity conditions 
insist that w = ±1) and the Schwarzschild-adS solution. To see this, suppose 
that w > everywhere. Then w cannot have a local minimum in the region 



< ti; < 1 nor a local maximum in the region w > 1 13 1. For regular 
solutions, w = 1 at the origin, and if it is initially increasing, it is necessary for 
w to first reach a local maximum if it is to cross w = 1 again, which cannot 
happen. Similarly, if w is initially decreasing, then, because w has no zeros, 
it cannot have a local minimum which is necessary if w is to cross 11! = ! 
again. For black hole solutions, the argument is dependent upon the value of 
w{rh). If w{rh) = 1, then w = 1 for all r, and the black hole geometry is 
Schwarzschild-adS. If w{rh) > 1, then initially w is increasing, and the same 
argument used above for the regular solutions applies. If w{rh) < 1, then 
initially w is decreasing, and once again w cannot subsequently increase to 
cross w = 1. Exactly the same argument works in both cases if w is negative 
everywhere rather than positive. 

In addition to the monopole solutions described in this section, dyonic 
soliton and black hole solutions also exist [Q] , again in contrast to the situation 
in asymptotically flat space, when the electric charge must be zero [|l2|. Due to 
the non-vanishing electric field, the stability analysis of these solutions is more 
complex (see, for example, p5[ | for the stability analysis of hairy black holes 
with non- vanishing electric field), and in this paper we shall focus only on the 
monopole solutions. However, we would anticipate that those dyonic solutions 
for which w has no nodes would also be stable in the odd-parity sector, at least 
when the electric field is small enough. 
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3 The pulsation equations 



We therefore consider linear fluctuations of a static and purely magnetic soliton 
or hairy black hole solution in Einstein- Yang-Mills theory with a cosmological 



constant. In a recent letter |11], we have shown that, at least in the pure EYM 
case, such fluctuations can be described by a symmetric wave equation for the 
linearized extrinsic curvature and the linearized electric field. In this section, 
we review the results obtained there, and generalize to include a cosmological 
constant. We discuss only the main steps in the derivation of a symmetric 



and hyperbolic formulation for a static background and refer to [16| for more 
details. 

One starts with the ADM equations of EYMA theory, where the metric 
and the gauge potential assume the form 

g = -a^dt^ + gij{dx' + (3'dt){dx^ + P^dt), 
A = -^adt + Ai{dx' + /3'dt). 

(<I> and Ai are both Lie algebra valued.) On the background, the slicing T,t is 
adapted to the staticity, i.e. the shift f3 and the time derivative of the 3-metric, 
Qij vanish. As a consequence, the extrinsic curvature tensor 

Kij = ^ {Sij - Lpgij) (2) 

vanishes on the background. Similarly, the electric YM field, defined in terms 
of the field strength F = dA + AAA and the future pointing normal unit vector 
field orthogonal to the slices, n = {dt — /?)/«, by 

Ei = Fif.n'', (3) 

is zero for static and purely magnetic solutions. Since Kij and vanish 
on the background, the tensors 6Kij and SEi, describing linear fluctuations, 
are invariant with respect to both infinitesimal diffeomorphisms within the 
slices and infinitesimal gauge transformation of the gauge potential. Hence, 
it is natural to look for a wave equation in terms of these "vector-invariant" 
quantities. In order to get such an equation, one then differentiates some 
evolution equations with respect to the time coordinate t and uses the linearized 
version of equations (§) and (|^) in order to eliminate time derivatives of 6gij 
and the gauge potential 6Ai . Next, one uses the linearized momentum and 
Gauss constraint equations and spatial derivatives thereof in order to make the 
spatial operator both elliptic and (formally) self-adjoint (the relevant scalar 
product is given below). Finally, one makes use of the freedom in choosing a 
(space-dependent) reparametrization of the time coordinate in order to impose 
the maximal slicing condition, 

6K\ = 0. 
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As a result, a hyperbolic and symmetric wave equation is obtained from the 
following combinations of the field equations, 

Af ^) = -a^tSiD^'Fi^) + ^Dia^SiDt'Fo^) + 2a^F^^Eko, 

where E^^^ = G^u - SiiGT^y = G^u - SttGT^u^^ + g^u^- Here and in the 
following, all quantities with a bar refer to the background 3-metric gij and 
the background magnetic potential Ai . Latin indices are raised and lowered 
with gij, while the index zero refers to the unit normal n. Finally, D = d+ [A\ . ] 
denotes the covariant derivative with respect to the gauge potential. 
The wave equation reads 

= Ai,- = ASr^^+Ajr*^ 

= (4) 

where Ajj denotes the trace-less part of A^. In terms of the vector-invariant 
quantities 

Lij = aSKij, A = Sa, £i = aSEi , 
one finds, after using the background equations E^i^ = and D>^F,yn = 0, 

A,^r^ = aLLij + 4V(, - 4a(,V'=L,)fe - 2aV'^ (^) L,)^ 



—V^iO^Vj) (^^^ + ^gij (^-^V\aa^)Lki + AA- R^oA^ , 



A^"""*^ r 1 1 



Af"^) = nESi + 4GaTi(P^£i)Fij 



+2aP''VkLij - '^Lkj&{a^F^') - aF^V^ J^^) , 

where Tr stands for an Ad-invariant scalar product on the Lie algebra and 
where all second partial derivatives are contained in the hyperbolic operators 
□l and De, defined by 



a 



OLLij = -d/ - VaVfc Lij + 2aR%Lj)k - 2aRkiijL^\ 
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Note that the spatial parts of the operators defined in are symmetric with 
respect to the scalar product 



((L«,£:«),(L(2),£:(^))) 

y'a^'l^I^ + 2GTr {g^^8\'hf]\ V^d'x. (5) 



s 



The constraint equations are the linearized momentum constraint, 



= a6Eio = aV^ (^^^ - 2GTr (F/,) , (6) 



and the linearized Gauss constraint 



= -a6{D^'Fo^) = aD^ 



Additional constraints involving also perturbations of the metric and the gauge 
potential themselves are the Hamiltonian constraint and all evolution equa- 
tions, which we had differentiated in time in order to construct the wave oper- 
ator. 

Since we adopt the maximal slicing condition, there is an elliptic equation 
for A, which is obtained from the trace of the tensor Ajj , 

(A - i^oo) A = [a'Lki + GTT{aFik8')] - 2GaTr [f^^'P^Lu] , (7) 

where the momentum constraint equation and the background equations 
have been used in order to simplify the equation. Here Rqq is the 00-component 
of the Ricci tensor. From the background equations, one finds that Rqq = 
GTi{F^^ Fj^i) /2 — A, and therefore, for a negative A, the operator on the left- 
hand side of equation (|^) is negative and the equation is solvable. As we have 



argued in |16|, this equation decouples when one projects the wave operator 
onto the momentum constraint manifold. In the case of odd-parity perturba- 
tions on a spherical symmetric background, as we are going to study below, 
equation is in fact void since it is a scalar equation. 



4 Odd-parity fluctuations 

4.1 The form of the perturbation equations 

We now apply the curvature-based perturbation theory in order to study the 
linear stability of the solutions discussed in section ^. Since the background 
is spherically symmetric in this case, it is convenient to expand the linearized 
extrinsic curvature and electric YM field in terms of spherical tensor harmonics. 
Perturbations belonging to different choices of the angular momentum numbers 
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a. and m decouple. Furthermore, the tensor harmonics can be divided into 
parities. Here, we consider only odd-parity (axial) perturbations. The even- 
parity (polar) case will be subject to a future article. Since the perturbations 
of all scalar quantities vanish for odd-parity perturbations, the variation of 
the lapse and the trace of the extrinsic curvature do not appear and therefore 
the amplitudes parametrizing 5Kij and 5Ei are fully coordinate- and gauge- 
invariant. 

How to expand Lij and £i in spherical tensor harmonics and how to find 
the corresponding expansion of the wave operator is explained in Appendix 
We introduce a coordinate p such that dp = NSdr ■ Since we are dealing 
with geometries which are asymptotically adS, p does not tend to infinity as 
r ^ oo, but instead p tends to a finite value pmax- We shall use the notation po 
to denote the lowest value of p, which is at the origin for solitonic solutions, 
and — oo at the event horizon for black holes. Using this "tortoise" coordinate, 
the resulting wave equation is 



where U = {h, a, b, k, c, d, e), with h and k parametrizing metric perturbations 
and a, b, c, d and e parametrizing the YM field. The matrices A and S are 
antisymmetric and symmetric, respectively, and have the block structure 




(8) 




where Ai and Si are 3x3 matrices, whereas A2 and S2 are 4x4 matrices. 
Defining 




the non-vanishing matrix elements of Ai and A2 can be written as 




The matrix S is given by 



sym. 



sym. 



Si = 7' 



p V 



+ [A + 2/] 
—2pw 



sym. 




7^ / 




sym. 



sym. 




V\v [A 





sym. 



sym. 




if^' -f] J 



9 



7' 



\ 



2VXlp 
p 

V2GXju 



\ 

-2(^10) p — 2\/TJ'juv 

-\/2//7p 
V2Xjp ) 



where /i^ = ^(^+ 1) and A = /i^ — 2 depend on the angular momentum number 
i. The constraint equations are found to be 









dp + 2VGu b - VXjk 



JVC, 



fi-y 



c — V2jw d, 



7dp{^b^ +JWC + -^7 d - y ^ 7 e. 



(9) 



How to solve the initial value problem in terms of gauge-invariant quantities 
is discussed in Appendix^. A nice fact is that the cosmological constant does 
not appear explicitly in the above system, it appears only via the background 
quantities S, N and w. As a consequence, the stability can be discussed on 
the same lines as for asymptotically flat solutions Q. Below, we study per- 
turbations which are smooth and which vanish at the boundary points p = po 
(r = or r = r/i as applicable) and p = Pmax ir = oo). The spatial part of the 
wave operator defined in (^) is symmetric for such perturbations, and since the 
operator is real, self-adjoint extensions exist. The spectrum of the self-adjoint 
extension which corresponds to physical boundary conditions decides the lin- 
ear stability of the background solution. It is, however, beyond the scope of 
this article to determine the self-adjoint extensions. We will rather base our 
stability argument on an energy estimate. 



4.2 The projection onto the constraint manifold and stability 

In order to study the linear stability of the background solutions, one has to 
consider the time evolution of perturbations C/, which are restricted to the 
constraint manifold defined by (|9|). Our next aim is therefore to separate 
the purely dynamical degrees of freedom from the constraint violating modes. 



Generalizing ideas introduced by Anderson et al. |17], we define the constraint 
variables Uc by the right-hand side of (^). We are then looking for a dynamical 
variable, Up , say, such that the wave equation assumes the form 



d 



dp + 



Ur 



0, 



(10) 



when expressed in the new variables {uc,Up). The operators Vc , Vj 
are permitted to have up to first order spatial derivatives only. The components 



pc and Vp 
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in Vcp have to vanish in order for the constraint variables to evolve. On the 
constraint manifold, Uc = 0, we then get 

- + Vp] Up = 0, 

which is a wave equation for the dynamical variables Up . Of course, it would 
be nice if the transformation can be found in such a way that the new potential 
Vp is automatically symmetric. In our case, it turns out that this is possible. 

As a simple example we first decouple the dynamical modes from the con- 
straint violating ones in the vacuum case where the YM amplitudes a, 6, 
c, d and e vanish, and where w = 1. The idea is to find a first order linear 
transformation of the form 

such that the spatial part of the wave operator factorizes. 




B^B, (11) 



where B^ = —dp + C"^ is the (formal) adjoint of B. The desired wave equation 
( [l0| ) is then obtained upon applying the operator B to the left of the original 
wave equation (|8|). This yields 



d? + BB^ 



0. 



Since BB^ is symmetric, the fact that Vcp must vanish implies that Vpc has to 
vanish as well. Provided that the factorization (|ll|) can be found, the constraint 
and dynamical variables can therefore be decoupled completely. 

It is far from obvious that such a transformation exists, since equation (|Tl| ) 
automatically implies that the spatial operator is positive. Nevertheless, it 
turns out that the factorization can be found in our case. Indeed, by making 
the ansatz 



-^7 



V -VXi A J 

where the first row of the matrix has been chosen such that Uc is defined by the 
right-hand side of the first equation in (P) and the second row such that there 
are no first order derivatives in B^B, equation (|ll]) is fulfilled if ^ = —-y- The 
transformed spatial operator is 

BB^ = -d''+ I VT/pp ^ 



p 



r{^)pp + l'^ 
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Using r{l/r)pp = 7^(2 — Gm/r) for a Schwarzschild background, we recognize 



that the dynamical variables are governed by the Regge- Wheeler equation 1 15 ] . 
Since the spatial operator has the form B^B, it is positive on the space of 
smooth perturbations with compact support. Therefore, in the odd-parity 
sector, the linear stability of the Schwarzschild black hole can be established 
without needing the explicit form of the Regge- Wheeler potential. This sort of 
topological argument will be important below, since the background solutions 
in EYM theory are not known in closed form. 

We now show how to decouple the constraint and dynamical variables for 
the full EYM wave equation d^. Here, the ansatz is 



Uc 



BU, B = dp-A + Cs, 



where the antisymmetric matrix A is given above and is introduced in order to 
reproduce the first order derivatives in the wave operator in (P). The matrix 
Cs is assumed to be symmetric. Writing this matrix in block form. 



Ct 



Co 



we must have 



7 



7 






/ 



V 






2VGu 


\ 


_> 

7 









_> 

7 




— V 









-V2w 




w 


M 

V2 





in order for the constraint variables to be defined as the right-hand side of (^). 
Hence, only the symmetric 4x4 matrix C2 has to be matched. 

Equating B^B with the spatial part of the wave operator, —d^ + Adp + 
dpA + S yields the following equations. First, we find the equation 

-dpCi + cl + cfCt + 2A1C1 = Si + dpAi , 

which is a consistency condition that can be shown to hold automatically. Next, 
we have a linear algebraic equation for the unknown, symmetric matrix C2 , 



CfC2 — Q^, 

where Q = St + dpCt — CjCi — A2Ct ■ Explicitly, we find 



(12) 





V GX U —Wn 



^fiu 



V2 
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Finally, we get a differential equation for C2: 

-dpC2 + Cl + [A2,C2\ 

where T = S2 - CtCl + Al , so that 



(13) 



V 



^ + 2Gu^ 



y/2Gun 



sym. 



sym. 
sym. 



V2 



7 7- 



,2rA 



7' 



sym. 
sym. 
sym. 



2Gu^ J 



Equation (|13| ) can be shown to be consistent with the linear equation (|l 
More precisely, if we set C = C2Cf — Q, we can show that the differential 
equation for C2 implies that 

-dpC + (C2 + A2)C - C{2Ai + Ci) = 0. 

Suppose that we found a solution to the equations (12) and (13). Then, 
the wave equation (^) is equivalent to the wave equation 

+ sst) V = (5^2 - + Adp + dpA + s')v = o, 

where V = BU . The symmetric potential S is given in block form by 

51 = Si + 2[CuAi] + 2dp{Ci + Ai), 
St = St + 2{CtAi-A2Ct) + 2dpCt, 

52 = S2 + 2[C2,A2]+2dpC2. 

Since Si and St do not depend on C2 , they can be computed directly. One 
obtains 



Si 

1' 



4Gu^ 



PP 



7 V7 



sym. 
+ A + 2/ 



V 



PP 



-2fj,w 



- - 

7 V7 



sym. 
sym. 
+ A + / 

PP 



\ 



and St = 0. The fact that St vanishes shows that - provided that we can 
solve the equation (13) - the constraint variable and the dynamical variables 
decouple from each other and are both governed by a symmetric wave equation. 
The evolution equation for the constraint variables is given by the symmetric 
wave equation 



dt' - dp' + Aidp + dpAi + Si]u, 



0, 



and therefore, initial data which satisfies Uc 
equations for all later times. 



Un 



will satisfy the constraint 
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If we have been able to factorize the spatial operator, stability of our solu- 
tions is then automatic. More precisely, one can show that for smooth pertur- 
bations with compact support, the energy expression 

pmax 

E = \ j {{U,U) + {BU,BU)]dp, 

PO 

is constant in time. Our boundary conditions are that all perturbations vanish 
at the origin (or event horizon) and at infinity. However, at this point, one can 
also require less restrictive boundary conditions. At the horizon, for example, 
it is sufficient to require that U and U are finite. At infinity {p = Pmax), one can 
impose the outgoing wave condition U = —BU, which means that there is no 
radiation coming from infinity. This condition implies that the energy cannot 
increase in time. Since E > 0, the kinetic energy must therefore be bounded 
and exponentially growing modes cannot exist. The remarks in Appendix ^ 
then show that the metric and the gauge potential themselves cannot grow 
exponentially. 



4.3 Factorization of the wave operator 

The result of the previous section is that we have shown the stability of our 
solutions once we have found a factorization of the spatial operator, that is. 



a solution of the equations (12) and (13). This is the subject of the present 
section. 

For the moment, we exclude the Schwarzscild-adS and the RN-adS back- 
grounds from the analysis below. These cases have to be treated separately. 
We start with the distinguished cases £ = and £ = 1. For i = 0, only the 
YM amplitudes a and d are present. The linear equation (jl^) then reduces to 
— \^wC2~^^ = V^Wp, with the unique solution C'^~^^ = —Wp/w. Equation 
([I3D is also fulfilled by this choice. Since we are interested in the fundamental 



solutions where w has no zeros, C2 '^^ is regular, and the transformed potential 
turns out to be 

S2 = 7 w 1 + 2 

V w 

which is the same as the one obtained in Q in the sphaleronic sector. The 
potential ^2 is positive, and has the asymptotic behaviour 

st'' ^ ^ as r 0, 
O2 — > U as r — > r/j, 
^2 — >--^(l + ^^oo) as r — > 00, 

the behaviour at r = being the same as for p-waves. Since 5^2^ is positive, 
linear stability follows, as discussed at the end of the previous subsection. 
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[Note that for solutions where w has zeros, Cg^""^ is singular. In this case, 
the method presented in |]l9| can be generalized and the existence of exactly n 
unstable modes established (n being the number of nodes of the function w).] 
For i = 1 the amplitudes k and e are absent. The corresponding rows 
and columns in the matrices Ct, C2, St, S2 and A2 have therefore to be 
removed. [The consistency conditions checked so far remain valid, since for 
i = 1, the entries in the removed rows and columns decouple from the entries 
in the other rows and columns.] The linear equation ([l^ ) turns out to have the 
unique solution 

^(i=l) ^ Wp / w -1 

^ 1 - w"^ \ -I w 

which also satisfies the differential equation (|l^). Thus, the constraint and 
dynamical variables decouple, and the dynamical variables are governed by the 
following symmetric wave equation: 

- + ~st'^) u, = 0, 

for 



where the symmetric potential is given by 




~At=i) _ 2! 2 + 4G^^^^ 4«; V 2u;2 ( l + w^ -2w 



2 ' y 4w 1 + J {l-w^)^\ -2w l + u-^ 

For solutions in which w has no zeros, the supersymmetric transformation is 
regular provided that 7^ ±1. However, as discussed in section ^, for solitonic 
solutions w = ±1 only at the origin, and for black holes, w"^ is never unity 
unless w = ±1, in which case the geometry is Schwarzschild-adS. The new 
potential is again positive and has a similar asymptotic behaviour as that for 
i = 0. Therefore, the stability follows also in this case. 

For solutions in which w does have zeros, the transformation is again regular 
provided that w 7^ ±1. In asymptotically flat space, w does not cross ±1 for all 
regular solitonic and black hole solutions, except that w = ±1 at the origin for 



solitons 1 12]. This is also the case for the majority of solitons and black holes 
in adS, so we can conclude that these solutions also have no unstable modes 
for i = 1. However, there are some regular solutions for which w has zeros and 
does cross ±1 away from the origin or event horizon for which we are not 
able to draw conclusions about the existence of unstable £ = 1 modes. 

We now turn to the generic case i > 2: The general solution to the linear 
equation (^) can be written into the form 

C2 = DiXo + T{p)X^)D, 
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where D = diag(n2, w.3, m4, ms) is a positive, constant matrix (the values of 
722, ... , m5, which are not relevant for what follows, are given in Appendix^). 



The symmetric matrices Xq and Xi are 



Xr 



X, 



\ 



-2Xfi'^^ + aVG/iwuv sym. 

—2\/2G fiw u 2wwp 

2V2Gfiu -2wp 

4V2G w^{l - w'^)u 2w'^Wp 

Aw'^v'^ sym. sym. sym. \ 
-2\j2GuP'v 2u? sym. sym 
2\j2Gvjv —2w 2 sym 

V2f2WV -f2W /2 



sym. 
sym. 


-2WWn 



sym. \ 

sym. 

sym. 

-f2WWp I 



1 f2 
2J1 



where /i = A + 2w^ and /2 = — 2?/;^. So far, the function T{^p) is arbitrary. 



Introducing this into equation ( |13D yields the following non-linear first order 
differential equation for T: 

- doT + AT"^ + BT + C 



with 



^2 AC 



0, 



(14) 



2G 



r 
2G 



(A + 2w^) + 1 



{w^ — l)wWp , 



(A + 2w^) + 2A + 



A 



(15) 



To summarize, it is sufficient to show that the single differential equation (jl^) 
admits a global solution with appropriate boundary conditions in order to show 
the stability of the evolution system (^), which is a wave equation for seven 
amplitudes. In the next subsection we shall complete the stability proof by 
showing the existence of a globally regular solution to (|l^. 

Finally, we turn to the stability of the Schwarzschild-adS and RN-adS black 
holes. The stability of the Schwarzschild-adS metric can be established on the 
same lines as above, but one has to take into account that for £ = 1, the 
matrix C2 is no longer uniquely specified by the linear equation (|l^), and 
one obtains a differential equation similar to (14). Alternatively, one can also 



use the perturbation formalism presented in Ref. |20|, since in the odd-parity 
sector, the cosmological constant does not appear explicitly. 

While the stability proof goes through for £ > 1, the RN-adS solution turns 
out to be unstable with respect to odd-parity radial perturbations of the non- 
Abelian part of the YM field which is parametrized by the amplitude d. When 
ti; = 0, the amplitudes a and d decouple, and d is governed by the equation 



0. 
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An easy way to find the instability is to use a trial function ||2T[. We define a 
sequence of functions d{p) as follows: 



d{p) 



for J = 1, 2, . . .. The function Z{p) G [0, 1] is defined by |21 



Z{p) 



and 



dZ ^ 
0< — <Q 
dp 



for p < -P - 1, 

1 for p > -P, 

for pG 1,-P], 



where P and Q are arbitrary positive numbers. The expectation value of the 
spatial operator yields, after some calculation, 



< 



91 
J 



2 , 
Ndr - :^\dr 



1 



while {d\d) < oo. This expectation value is negative provided we choose J 
to be sufficiently large. Therefore we conclude that there are exponentially 
growing modes and the RN-adS black hole is unstable. 



4.4 Global solutions to equation (^) 

It remains at this stage to show that equation (|T^ ) has globally regular solu- 
tions. We shall begin by considering only those solitons and black holes for 
which w has no zeros, since these are the solutions of most interest to us. 

In order to see that equation ( |T4[ ) admits global solutions, one rewrites it 
as a linear, second order differential equation using the transformation 



A z ■ 

Going back to the radial coordinate r, this yields 



+ 



1 d{NS) 
JlS dr 



A dr 



where 



p^XC 



^{X + 2w') + l 



[w 



^{X + 2w'^) + 2X + W 



dz ~ 

— + ACz = 0, 

or 



l)wWr, 

p^X fx 



(16) 



(17) 



Nr"^ \ 2 



— \- w' 



(18) 
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The equation (^) is regular everywhere except at the origin, infinity, and (if one 
exists) the event horizon r = rh, where there are regular singular points. Using 
the asymptotic expansions and the standard Frobenius method, the solutions 
for z near each of these points can be found. 

Near the origin, for globally regular (solitonic) solutions, the functions in 
([TtI) behave as: 

fi^XA = 2X^ + AX + 0{r^), 

so that the origin is a regular singular point, and the two linearly independent 
solutions have the form 

£+1 -£ 
z ~ , r . 

If, instead, we are considering black hole solutions with an event horizon at 
r = rh, then A and B are regular functions at the event horizon, while C 
diverges as (r — rh)^^. Therefore, again there is a regular singular point, this 
time the indicial equation for z has a repeated root of zero, so that the linearly 
independent solutions have the form: 

z = 0(l), 0{log{r-rh)). 

As expected, the negative cosmological constant means that the analysis at 
infinity is different to that in the asymptotically flat case |^]. Letting s = 1/r, 
the behaviour of the functions in (ITtI) is: 



NS = -A + 0(1/,), 

C = C^s^ + 0{s'), 

which implies that the point s = is, in fact, a regular point. As a consequence, 
the asymptotic forms of the linearly independent solutions for z at infinity are: 

z = 0(l), 0(l/r). 



The proof of the existence of global solutions to equation (17) can now pro- 
ceed as in the asymptotically flat case |^]. Consider firstly black hole solutions. 
We start near the event horizon with the solution having z = 0(1) (i.e. with 
no logarithmic term), so that z is positive for sufficiently small r — r^. At a 
maximum of z, where dz/dr = 0, from ( [l7| ) we have 



= -ACz. 

dr'^ 



18 



It is clear from ([T^) that „4 > always, so we need to investigate the sign of 
C. In the asymptotically flat case it was checked numerically that C < 0, 
but in our situation, we can prove this analytically for sufficiently large |A|. 
First fix £ = 2. It was proved in |8| (see also section ^ that, as |A| ^ oo. 



iV~0(|A|), ^-o(\A\ 



for all r. Hence, for sufficiently large |A|, both terms on the right-hand-side of 
(18) are vanishing, but the first term is o(|A|~^) and the second term 0(|A|~^), 



so the first term is vanishing more rapidly. Therefore, for sufficiently large |A|, 
we have that C < for all r, for this particular value of £. If we now replace 
i ^ then iJ?XC — > fi^XC + where a lengthy calculation (using the fact 

that C < for the original value of £) yields 



i)H£ + 2) 



-^{21^ + 6i) - 2X{2f + 6i) 



-8w\£ -!){£ + 2) i^ + l 



< 0, 



where A and /i are calculated with the original value of i. Therefore, as i 
increases C decreases, and is therefore negative for all values of i, for sufficiently 
large |A|. 

We conclude that z cannot have a maximum if it is positive, nor minimum 
if it is negative. Therefore, since z is positive close to the event horizon, it will 
remain strictly positive, and take the form 

z = Zoo + 0{r-^) 

at infinity, where z^o > 0. Since A is positive, and z has no zeros, using the 
definition ( [T6| ) we see that T is regular and exists globally with the asymptotic 
behaviour 

T ~ (r - rh) 



near the event horizon, and 



0(1) 



at infinity. As an example, the function T/S for a particular black hole solution 
with A = —100 is shown in figure |l|. 

For globally regular solutions, a similar argument holds, using the solution 
z ~ r^"*"^ near the origin. In this case the function T behaves like for small 
r but the behaviour at infinity is unchanged. 

A careful analysis reveals that again, the behaviour of the transformed 
potential 5^2 is similar to the £ = case: ^2 = 0(l/r^),0(r — r/i),0(l) near 
r = O,r/j,oo, respectively. 
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.12 I 1 ' ' 1 1 1 1 1 

012345678 

log(r) 

Figure 1: The function y = T/S for i = 2,3,4,5. The background solution is 
a black hole with A = —100, rh = 1, and w{rh) = 0.9. 



We have now completed the proof that those black holes and solitons which 
have no unstable modes in the odd-parity sector for spherically symmetric 
perturbations also have no unstable modes in the odd-parity sector when we 
consider non-spherically symmetric perturbations. In other words, we have 
proved the (odd-parity) stability of those solutions in which w has no zeros 
and |A| is sufficiently large. 

A similar analysis shows that the Bartnik-McKinnon solitons and the cor- 
responding black holes with hair have no unstable modes with £ > and odd 
parity Q. The only difference there, is that the background quantities imply 
that the function T behaves like 1/r near infinity. Instead of approaching a 
constant value, the transformed potential falls off as when r — > oo. 

The question, therefore, is whether our analysis can be extended to other 
values of |A|. The answer, surprisingly, is no. The crucial part of the work is 
the sign of the function C. For large |A| we have been able to show analytically 
that C is always negative. When A = 0, numerical analysis showed this to also 
be the case However, for small, negative A, we have found numerically that 
C can be positive and the solutions to ([T^ ) blows up (for example, for black hole 
solutions with A = —0.001, r/j = 1 and w{rh) = 0.642). Therefore our analysis 
does not extend and it may be that these solutions have additional instabilities 
(they are already unstable for spherically symmetric perturbations). In order to 
find such instabilities, one could use the nodal theorem for coupled, symmetric 
systems of Schrodinger equations proved in Ref. to count the number of 
bound states of the spatial operator. 
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5 Conclusions 



In this article we have studied the Unear stabihty of sohton and hairy black 
hole solutions of su(2) Einstein- Yang-Mills theory with a negative cosmological 
constant. We have applied a recently developed perturbation formalism which 
is based on curvature quantities. The advantage of this approach is that one can 
work with gauge-invariant quantities, which are governed by a wave equation 
whose spatial part is symmetric. As a consequence, we were able to decouple 
the constraint variables by a supersymmetric-like transformation, and to show 
analytically the linear stability of soliton and black holes with respect to odd- 
parity fluctuations. We stress that such a proof is unlikely to exist in a gauge- 
invariant metric formulation, since when non-Abelian gauge fields are coupled 
to the metric, the metric approach fails to yield a symmetric wave operator in 
a natural way. 

Our main result concerns those solutions in which the function w which 
determines the gauge field potential has no zeros and the cosmological con- 
stant is large and negative. These solutions are of particular interest as it 
was already known that they have no modes of instability under spherically 
symmetric perturbations. We have proved that this holds also for general, lin- 
ear, perturbations in the odd-parity sector. This result is significant because 
it shows that these static configurations possess no topological instabilities. 
Therefore the presence of a sufficiently large (negative) cosmological constant 
(which means that there is a large gravitational potential away from the black 
hole horizon, or the centre of the solitons) stabilizes the situation. This may be 
understood heuristically by analogy with the situation in quantum field theory 
in curved space. A black hole can be in thermal equilibrium with a bath of ra- 
diation at the Hawking temperature, however in asymptotically flat space this 
equilibrium is unstable. Stability can be restored either by placing the whole 
system in a box, or if the black hole instead resides in asymptotically anti-de 
Sitter space [^, provided the cosmological constant is sufficiently large and 
negative. In some sense, then, the structure of asymptotically adS space plays 
a similar role in our situation to placing the black hole or soliton in a box, 
which means that the gauge field cannot escape to infinity, as happens when 
the corresponding configurations in asymptotically flat space are perturbed 
[|4|. Further, b ecause of the above stability of the Hartle-Hawking state, the 
black holes we have studied in this paper are precisely those which are most of 
interest from the quantum field theory point of view. Therefore it is important 
to establish their classical stability before studying the properties of quantum 
fields propagating on these geometries. This will be the subject of future work. 

Classically, some open questions remain. Firstly, we need to investigate the 
stability of these solutions under non-spherically symmetric even-parity per- 
turbations. The even-parity is generally less amenable to analysis, since its 
properties depend crucially on the detailed structure of the equilibrium con- 
figurations. However, in the situation in which we are interested, when the 
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cosmological constant is large and negative, it was found in that, for spher- 
ically symmetric perturbations, the even-parity sector simplified considerably 
and stability was proved analytically in this case. Therefore it seems reason- 
able to suppose that simplification may be possible also for non-spherically 
symmetric perturbations. 

Next, we have not examined the zero modes of the pulsation equations, 
i.e. the stationary solutions to the perturbation equations. Similarly to the 
asymptotically fiat case we except to find slowly rotating and charged 
solitons and black holes, the latter ones generalizing the Kerr-adS metric. 

We shall return to these questions in a subsequent publication. 
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A Harmonic decomposition of the wave operator 

With respect to the background 3-metric 

g = h r^dO,\ 

odd-parity perturbations of the extrinsic curvature can be expanded according 
to 

Lrr = 0, LrB = Sb, LaB = rk 2V(^5b) . (19) 



Here and in the following, capital indices refer to coordinates on the 2-sphere. 
In the above equation, 5^ = are the transverse vector harmonics, and 
can be expressed in terms of the standard spherical harmonics Y = Y^"^ as 
Sb = {*dY)B, where a hat refers to the standard metric on the 2-sphere. 

Similarly, the electric YM field is expanded into su(2)-valued vector har- 
monics with odd parity (see [^] for details), 

£r = — 1= Xi H — X2 , 

rVN rVN 



Sa = cTrYA + dYTA + e(vAX2 + ^fi^YTAy (20) 

where in terms of the Pauli matrices a = (u*), = •£/ (2i), ta = •£/ (2i). 
Here Xi and X2 are the su(2)-valued harmonics 

X^=YTr, X2=g^''VAYTB. 



22 



The tensor harmonics in the expansions ([T^) and are chosen to be orthog- 
onal. After the rescahng 

h = nih, k = n2k, a = mia, b = mib, c = niic, d = mid, e = mie, 

where in terms oi iJp = £{£ + 1) and \ = fj? — 2, the coefficients are given by 

1 1 

ni = , — - , n2 



2/i2 

1 1 
mi = , 7712 = ms = , mi = , m^ 



the expansion is normahzed such that 

((L, 8) ,{L,£)) = I {h^ + + + 62 + c2 + ^2 + rf^ 



N 

where (., .) denotes the scalar product defined by (^). Note that for i = 1, 
^{aSb) ^-iid VaX2 + ^fi'^YTA vanish, and therefore, the amplitudes k and e 
do not exist in those cases. For i = 0, the function Y is constant and only the 
amplitudes a and d are present. 

An efficient way to perform the harmonic decomposition of the wave oper- 
ator is to compute the energy expression which corresponds to Ajj and A^^^^^, 

E = Egrav + EyM + Eint, 

where 



-2AV^Lij + 4GTr <| F^F^Lu + -PkiP'^' UjlA^ a^dx\ 



Eym = G J Tr + 2£)[^f^l • + & 

E^nt = -4G jTr!^U^F,^Dk£j + ^V^£kDi(^a^Fj'' 



For a spherically symmetric background with no electric field, non-vanishing 
background quantities are given by 

^ ArB = -h^rQAB. R%AB = 2(1 " ^)^\Amc. 
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Rrr = -^, RaB=[1-N -'^Nr)gAB, 

FrB = -Wr £\ta, FaB = {w^ " l)7"r EaB- 

Using this, the background equations and the expansions ( p^ ) and ([20|), the 
energy expression yields, after integrating over the spherical variables, 

E = \j [{U,U) + {dpU,dpU) + 2{U,AdpU) + {U,SU)yp, 

where U = {h,a,b,k,c,d,e), and where A and S are given in section ^. The 
radial coordinate p is defined by dp = SNdr, and (. , .) denotes the standard 
scalar product. The wave equation (|8|) follows directly from this energy func- 
tional. 



B The initial value formulation 



In this appendix, we give a gauge-invariant formulation of the initial value 
problem. In order to solve the linearized EYM equations, we have to take into 
account the Hamiltonian constraint and all evolution equations which we had 
differentiated in time. These equations cannot be described in terms of the 
linearized extrinsic curvature and the electric YM field alone, since perturba- 
tions of the 3-metric and the magnetic gauge potential appear (and not only 
their time derivatives). However, if the background is spherically symmetric, 
one can give a formulation in terms of gauge-invariant quantities which we 
introduced in an earlier article |20|. 

It turns out that the "missing" constraint equations (i.e. the linearized 
Hamiltonian constraint and the relevant evolution equations we had differenti- 
ated in time) are equivalent to the p-components of the equations (41), (42) and 
(43), and to equation (46) of Ref. ||2^ [The presence of the cosmological con- 
stant does not modify those perturbation equations.] Using the gauge-invariant 
amplitudes H, A, B and C defined there, these equations become 



2dth 



did 



dtb 



r 



4G7 



(w^ - l){Ap - wBp) 



WpC 



(1 



w 



2\2 



-H, 



7 



{p^ + 2w'^)Ap - 2p^wBp - 2wCp + 2wpC + fi^{l - w^)^ 



-2wAp + (^2 + y;2 _ i)Bp +Cp- w{l - w^)^ 



dt 



d+-p 



2 ;^ 



wdpAp + 2wpAp + 72[/i2 ^yj-i -i]c 



Cpp + p?"w^ 



P J.2 



(21) 



Furthermore, one has the relation between the gauge-invariant amplitudes of 
Ref. [|Q] and the curvature-based amplitudes /i, k, d, ... introduced in appendix 
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A of this article. For ^ > 2, one has 
Ht = -2rk, 

Bt = e, 
C = -d- \n'^e + wAt. 

The initial value problem for I >2 can be solved as follows. First, we choose 
any functions Hp = H^p \ Ap = A^p \ Bp = B^^^ and C = C^^^ on an initial time 
slice, St=o, say. Next, one solves the momentum constraint equations for 
U = {h, o, b, k, c, d, e) on the initial time slice. A convenient way to do this is to 
freely specify the functions c, d, e and k and to compute h, a and b using (^). 
Then, the time derivative of U on Eg is consistently given by equations (pl|) and 
the time derivatives of the momentum constraint equations. The amplitudes 
U are then evolved using the symmetric wave equation (^). Finally, the gauge- 
invariant amplitudes H, A, B and C, parametrizing the metric and the gauge 
potential are obtained from (^) after integration over t, where the integration 
"constants" are given by Hp'\ Ap'\ Bp'^ and C^^\ For i = 1, the perturbation 
equations can be solved in a similar manner, but one has to take into account 
that C = and that the relation (^) is different in that case. 

Since we have shown that the pulsation equations admit no solutions which 
grow exponentially in time when |A| is large enough, the relations (^) show 
that the same must hold for the gauge-invariant quantities H, A, B, and C, 
parametrizing the metric and the gauge-potential. This completes our stability 
analysis. 
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